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ABSTRACT 


The  purpose  of  this  technical  note  is  to  present  a  new 
algorithm  to  obtain  real-time  information  for  electronic  support 
measures  for  wideband,  low  probability  of  intercept  radar  and 
communication  signals.  This  algorithm  is  Cyclic  Cross 
Correlation (CCC) .  It  allows  the  detection,  separation  and 
characterization  of  spread  spectrum  signals  having  overlapped 
spectra,  the  same  carrier  frequency  and  the  same  time  of  arrival. 

An  optoelectronic  real-time  implementation  of  the  CCC  for  large 
bandwidth,  spread  spectrum  signals  was  constructed  at  the  Defence 
Research  Establishment  Ottawa  using  an  optical  Time-Integrating 
Correlator (TIC) .  A  description  of  the  real-time  optical 
implementation  of  the  CCC  with  a  TIC  and  the  associated 
experimental  results  are  presented.  An  assessment  of  the  potential 
of  this  technique  for  application  to  radar  and  communication  ESM  is 
also  included. 


RESUME 

Cette  note  technique  prSsente  un  nouvel  algorithme  congu  pour 
extraire,  en  temps  rSel,  de  1' information  de  signaux  de  radar  et  de 
communication  &  large  bande  passante  et  5  faible  probability 
d' interception  pour  les  mesures  de  soutien  Slectronique.  II  s'agit 
de  1'intercorrSlation  cyclique  (IC) .  Cet  algorithme  permet  de 
separer  et  de  caracteriser  des  signaux  a  spectre  Stales  ayant  la 
mSrae  porteuse,  le  mSme  temps  d'arrivSe  et  des  spectres  superposes. 
La  mise  en  oeuvre  optoSlectronique,  en  temps  reel  de  1'IC,  par  un 
Correlateur  A  Integration  Temporelle  (CIT)  a  ete  realisee  au  Centre 
de  recherche  pour  la  defense,  Ottawa,  pour  des  signaux  A  spectre 
etale  a  large  bande  passante.  Ces  travaux  ainsi  que  les  r§sultats 
experimentaux  obtenus  sont  presentes.  Une  evaluation  du  potentiel 
de  cette  technique  pour  application  aux  mesures  de  soutien 
electronique  pour  le  radar  et  les  communications  est  prSsentee. 


iii  . 


EXECUTIVE  SUMMARY 


Cyclostationary  processing  techniques  are  emerging  as  powerful 
analysis  tools  for  modern  radar  and  communications  spread  spectrum 
Low  Probability  of  Intercept  (LPI)  signals  that  use  pseudo  random 
sequences.  These  signals,  either  spread  with  sequences  only  or 
hybrids  using  time  or  frequency  hopping  combined  with  sequence 
spreading,  are  characterized  by  large  bandwidths  and  low  power 
spectral  densities.  This  makes  their  detection  by  linear  receivers 
such  as  spectral  analysers  very  difficult.  Non-linear  receivers 
using  cyclostationary  processing  have  been  proposed  as  a  better 
approach  to  the  detection  and  interception  of  these  spread  spectrum 
LPI  signals  because  of  their  ability  to  sort  the  signals  from  the 
noise  by  using  new  parameters.  However,  the  computer  power 
required  for  a  real  time  implementation  of  these  algorithms  for 
large  bandwidth  signals  is  prohibitive,  even  with  the  best  digital 
systems  available.  Fortunately,  it  was  recently  demonstrated  at 
NRL  that  acousto-optic  wideband  correlators  substantially 
outperformed  their  digital  counterparts.  So  harnessing  the  power 
of  an  optical  correlator  for  the  implementation  of  cyclo  stationary 
signal  processing  is  an  exciting  challenge  that  could  lead  to 
significant  advancement. 

The  purpose  of  this  technical  note  is  to  present  a  new 
algorithm  to  obtain  real-time  information  for  electronic  support 
measures  for  wideband,  low  probability  of  intercept  radar  and 
communications  signals.  This  algorithm  is  Cyclic  Cross 
Correlation (CCC) .  It  allows  the  detection,  separation  and 
characterization  of  spread  spectrum  signals  having  overlapped 
spectra,  the  same  carrier  frequency  and  the  same  time  of  arrival. 
A  real-time  optoelectronic  implementation  of  the  CCC  for  large 
bandwidth,  spread  spectrum  signals  was  performed  at  the  Defence 
Research  Establishment  Ottawa  using  an  optical  Time-Integrating 
Correlator (TIC) .  A  description  of  the  real-time  optical 
implementation  of  the  CCC  with  a  TIC,  and  the  associated 
experimental  results  are  presented  together  with  an  assessment  of 
the  potential  of  this  technique  for  radar  and  communication  ESM. 
The  definition  of  the  CCC  is  introduced  in  Section  2  and  followed 
by  a  description  of  a  TIC  in  Section  3.  A  description  of  the 
implementation  of  the  CCC  with  a  TIC  appears  in  Section  4  with 
experimental  results  in  Section  5. 

The  real-time  implementation  of  the  CCC  has  great  potential 
for  the  detection,  characterization  and  direction  finding  of  spread 
spectrum  signals  in  a  crowded,  noisy  environment  where  spectrum 
overlap  is  common  and  more  than  one  signal  may  be  associated  with 
a  particular  carrier  frequency  and  time  difference  of  arrival. 
These  features  could  be  incorporated  to  design  and  implement 
advanced  Electronic  Support  Measures  against  LPI  spread  spectrum 
signals  used  by  some  modern  radar  and  communications  systems. 
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1 . 0  INTRODUCTION 


Cyclostationary  processing  techniques  are  emerging  as 
powerful  analysis  tools  for  modern  radar  and  communications  spread 
spectrum  Low  Probability  of  Intercept  (LPI)  signals  that  use  pseudo 
random  sequences.  These  signals,  either  spread  only  with  sequences 
or  hybrids  using  time  or  frequency  hopping  combined  with  sequence 
spreading,  are  characterized  by  large  bandwidths  and  low  power 
spectral  densities.  This  makes  their  detection  by  linear  receivers 
such  as  spectral  analysers  very  diff icult[l] .  Non-linear  receivers 
using  cyclostationary  processing  have  been  proposed  as  a  better 
approach  to  the  detection  and  interception  of  these  spread  spectrum 
LPI  signals  because  of  their  ability  to  sort  the  signals  from  the 
noise  by  using  new  parameters.  However,  the  computer  power 
required  for  a  real  time  implementation  of  these  algorithms  for 
large  bandwidth  signals  is  prohibitive,  even  with  the  best  digital 
systems  available.  Fortunately,  it  was  recently  demonstrated [ 2 ]  at 
NRL  that  acousto-optic  wideband  correlators  can  substantially 
outperform  their  digital  counterparts.  The  particular 
optoelectronic  correlator  presented  in  the  NRL  publication  had  a 
processing  power  of  10-12  Gflops  and  was  equivalent  to  25  to  60  VME 
boards.  So  harnessing  the  power  of  an  optical  correlator  for  the 
implementation  of  cyclo  stationary  signal  processing  is  an  exciting 
prospect  that  could  lead  to  significant  advancement. 

The  purpose  of  this  technical  note  is  to  present  a  new 
algorithm  to  obtain  real-time  information  for  electronic  support 
measures  for  wideband,  low  probability  of  intercept  radar  and 
communications  signals.  This  algorithm  is  Cyclic  Cross 
Correlation (CCC) .  It  allows  the  detection,  separation  and 
characterization  of  spread  spectrum  signals  having  overlapped 
spectra,  the  same  carrier  frequency  and  the  same  difference  of  time 
of  arrival.  An  real-time  optoelectronic  implementation  of  the  CCC 
for  large  bandwidth,  spread  spectrum  signals  was  constructed  at  the 
Defence  Research  Establishment  Ottawa  using  an  optical 
Time-Integrating  Correlator  (TIC) .  Optical  TICs  are  analog  computers 
designed  primarily  to  perform  the  correlation  of  large  bandwidth, 
long  data  streams.  Applications  such  as  DNA  analysis [3], 
processing  for  code  division  multiple  access  signals  and  searches 
in  large  unstructured  data  bases  have  already  been  proposed  for 
TICs.  A  description  of  the  real-time  optical  implementation  of  the 
CCC  with  a  TIC  and  the  associated  experimental  results  are 
presented.  An  assessment  of  the  potential  of  this  technique  for 
radar  and  communications  ESM  applications  is  also  included.  It  is 
demonstrated  that  the  optoelectronic  implementation  of  the  CCC 
allows  the  real-time  calculation  of  the  CCC.  Consequently,  the 
separation  and  characterization  of  spread  spectrum  signals  with 
overlapping  spectra  that  share  the  same  carrier  frequency  and  the 
same  time  difference  of  arrival  is  achieved.  An  assessment  of  the 
potential  of  this  technique  for  ESM  on  spread  spectrum  radar  and 
communications  signals  is  also  included  in  the  conclusion. 
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2.0  INTRODUCTION  TO  CYCLOSTATIONARY  SIGNAL  PROCESSING 

Many  signal  processing  techniques  assume  that  both  the 
signal  and  the  noise  are  stationary.  However,  man-made  signals 
such  as  radar  and  communications  signals,  usually  contain 
parameters  that  vary  periodically  with  time.  The  assumption  of 
stationarity  is  then  not  valid.  In  order  to  better  deal  with  these 
non  stationary  signals,  W.  Gardner  et  al. [2,4-11]  have  developed 
the  theory  of  cyclostationary  signal  processing  in  which  the 
periodicities  of  the  signals  are  recognized  and  exploited  by 
non-linear  receivers.  Therefore,  the  non-linear  receivers  based  on 
cyclostationary  processing  are  better  suited  for  the  detection  and 
interception  of  non  stationary  signals. 

Cyclostationary  signal  processing  is  a  field  of  intense 
activity  and  numerous  publications  are  available,  including 
tutorials [9]  and  textbooks [  10 ] .  Applications  have  been  proposed  in 
the  field  of  interception,  signal  feature  detection  and  direction 
finding  of  weak  man-made  signals [4-13 ] .  It  also  has  been 
recognized [14-15]  that  many  specialized  transforms  (often  developed 
empirically)  performed  in  the  time  or  the  frequency  domain,  such  as 
the  ambiguity  function,  the  Wigner-Ville  or  the  Choi-Williams 
distributions  are  cyclostationary  processes.  Optical 
implementations  have  been  proposed  for  some  of  these  processes  such 
as  the  chip  rate  detector [16] ,  the  ambiguity  function [17-19]  and 
the  cyclic  spectrum[20] . 

2.1  The  Cyclic  Cross  Correlation (CCC) 

The  particular  cyclostationary  process  that  will  be 
focussed  on  here  is  the  CCC.  The  CCC  function,  Rsis2(T)'  defined 
by 


R 


a 

sls2 


(T) 


lim  1  rT' 2e 

“  r-oo  T  J_T/2si 


(fc  +  t/2)  s2  (  fc-t/2)  e~i2*at  dt 


(1) 


This  is  a  qrosscorrelation  function  where  a  is  the  cyclic  frequency 
and  where  s1(t+r/2)  and  s2(t-r/2)  are  the  signals  received  by  two 
different  antennas.  The  Fourier  transform  of  the  CCC,  the  cyclic 
spectral  density  function  S®ls2(f),  is  defined  by 

S*ls2(f)  =  f~  Rsis2( t)  dt  (2) 

J  -o# 


For  a=0  it  reduces  to  the  conventional  power  spectral  density 
function  S°ls2(f). 
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Equation  1  can  be  interpreted [21]  as  the  correlation 
between  the  RF  signal  x(t+r/2)  received  by  one  antenna  and  the 
frequency  shifted  version  x(t-r/2)  e-l2,rat  of  the  same  RF  signal 
received  by  another  antenna.  It  can  be  demonstrated  that  for 
cyclostationnary  signals,  a  correlation  exists  between  frequency 
shifted  versions  of  the  signal  when  the  frequency  shift  is  equal  to 
one  of  the  periodicities  of  the  signal.  This  results  from  the 
periodical  nature  of  the  signal  in  time[9].  This  property  is  an 
important  means  of  distinguishing  cyclostationary  from  stationary 
signals  and  noise.  It  also  allows  the  detection  of  the  signals  and 
the  characterization  of  their  parameters. 

2.2  Illustration  of  the  Cyclic  Cross  Correlation 

In  order  to  help  visualize  the  CCC  and  to  understand  its 
properties  in  terms  of  discriminating  signals,  a  computer  generated 
image  has  been  produced  (see  Figure  1)  that  shows  the  CCC  of  three 
signals  with  overlapping  spectra  which  share  the  same  Time 
Difference  of  Arrival  (TDOA)  and  the  same  carrier  frequency.  It 
illustrates  the  CCC  of  the  three  spread  spectrum  signals  which  have 
been  used  for  the  first  set  of  experimental  results  of  Section  5.1 
of  this  Technical  Note.  The  geometry  of  the  location  of  the 
emitters  of  the  three  signals  is  illustrated  in  Figure  2  and  their 
spectra  are  illustrated  in  Figure  3.  These  signals  include  two 
Binary  Phase  Shifted  Keyed  (BPSK)  signals  having  the  same  TDOA  of 
t= 0,  the  same  carrier  frequency  but  chip  rates  of  2  MHz  and  5  MHz 
respectively.  The  third  signal  is  a  30  MHz  broadband  source  of 
Additive  White  Gaussian  Noise  (AWGN)  located  at  a  TDOA  of  t=-1.46 
M s.  Figure  3a)  and  3b)  illustrate  respectively  the  spectra  of  the 
2  MHz  and  of  the  5MHz  BPSK  signals.  Figure  3c)  illustrates  the 
spectrum  of  the  noise  and  Figure  3d)  the  spectrum  of  the  sum  of 
the  three  signals.  The  spectra  of  the  three  signals  are  overlapped 
and  the  BPSK  signals  have  the  same  carrier  frequency.  The  large 
amount  of  noise  present  makes  the  assessment  of  the  situation  very 
difficult  if  only  conventional  spectral  analysis  of  the  added 
signals  is  considered. 

The  CCC  illustrated  in  Figure  1  shows  many  interesting 
features  that  are  listed  below.  The  horizontal  plane  of  the  CCC 
(Figure  1)  contains  the  axis  of  the  cyclic  frequency  a  and, 
perpendicular  to  it  is  the  TDOA  r.  The  amplitude  of  the  CCC  is 
represented  (with  a  linear  scale)  on  the  vertical  axis.  It  should 
be  noted  that: 

1)  the  plane  a=0  contains  the  conventional  cross 
correlation  of  the  signals.  The  peak  N  for  the  AWGN  at 
t=-1 . 46  jus  and  the  peak  A  from  the  two  BPSK  signals  at 
t=  0  are  separated  by  their  different  TDOA. 
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Figure  3:  Spectra  of  the  signals  of  the  three  emitters: 
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2)  it  is  impossible  to  separate,  in  the  plane  a=0,  the 
correlation  peaks  of  the  two  BPSK  signals  (peak  A)  as 
they  share  the  same  TDOA  of  0  jus. 

3)  the  correlation  peak  of  the  AWGN  consists  only  of  the 
thin  peak  N  a  r=- 1.46  jus  and  a=0.  The  small  width  of 
the  peak  is  a  direct  consequence  of  the  large  bandwidth 
of  the  AWGN. 

4)  in  the  plane  r=0,  four  peaks  appear.  Peak  A,  at  a=0, 
has  already  been  mentioned  and  corresponds  to  the 
overlap  of  the  autocorrelation  of  the  two  BPSK  signals 
that  share  the  same  TDOA. 

5)  also  at  t=0,  peak  B  at  a=2  MHz  corresponds  to  the  CCC 
of  the  2  MHz  BPSK  signal  and  the  peak  C  at  a=4  MHz  is 
a  harmonic  of  the  CCC  of  the  same  signal. 

6)  finally,  peak  D,  at  r=0  and  at  a=5MHz,  is  produced  by 
the  CCC  of  the  5  MHz  BPSK  signal. 

Figure  1  clearly  demonstrates  that  the  CCC  has  the 
capability  to  detect  and  separate,  in  the  presence  of  a 
considerable  amount  of  noise,  spread  spectrum  signals  having  the 
same  carrier  frequency,  overlapping  spectra  and  the  same  TDOA. 
These  features  are  of  prime  interest  in  the  design  of  ESM  systems 
for  LPI  radar  and  communications  signals. 


3.0  TIME-INTEGRATING  CORRELATORS (TIC) 

TICs  are  analog  optical  computers  designed  to  perform  the 
correlation  of  two  signals.  There  are  many  ways  to  build  optical 
time-integrating  correlators.  They  are  well  documented  in  the 
literature  [22-41]  and  good  review  papers  are  available[29, 35,41] . 
The  various  factors  affecting  the  performance  of  these  systems  have 
also  been  extensively  studied[42-5l] .  A  brief  review  of  the 
principle  of  operation  of  these  correlators  is  presented  here, 
based  on  the  tandem  architecture  described  in  [52]  (Figure  4  and 
5).  The  emphasis  is  on  the  characteristics  and  parameters  that 
have  an  impact  on  the  implementation  of  the  CCC. 

3.1  Description  of  the  Optical  Processor 

The  first  operation  performed  by  a  TIC  is  the 
transformation  of  the  RF  signals,  a(t)  and  b(t) ,  to  be  correlated 
into  light  modulated  beams.  For  the  implementation  of  the  CCC,  the 
transmitted  RF  signals  should  be  received  by  two  separate  antennas. 
The  RF  signals  from  each  antenna  are  individually  applied  to  one  of 
the  light  modulators  (Bragg  cells) .  They  consist  of  a 
piezoelectric  transducer  attached  to  a  crystal  in  which  acoustic 
waves  are  generated.  As  the  acoustic  waves  propagate  in  the  Bragg 
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cell  crystals  and  a  moving  grating  of  changing  indices  of 
refraction  is  formed.  The  gratings  in  the  Bragg  cells  are 
illuminated  by  two  expanded  laser  beams  originating  from  the  same 
laser.  Each  beam  interacts  with  the  acoustic  waves  of  only  one 
Bragg  cell.  Some  of  the  incident  light  is  diffracted  by  the 
gratings  and  therefore,  the  information  contained  in  the  RF  signals 
is  simultaneously  transferred  to  the  diffracted  laser  beams.  The 
signals  propagating  in  each  Bragg  cell  are  physically  oriented  in 
opposite  directions,  a  condition  needed  for  the  correlation.  The 
gratings  are  imaged,  with  a  lens,  onto  a  detector  array  where  they 
remain  counterpropagating.  The  5000-element  detector  array 
performs  a  time  integration  on  the  coherent  addition  of  the  two 
images  of  the  signals. 

If  the  variable  z  is  the  distance  along  the  Bragg  cells  and 
their  images,  and  if  the  origin  z=0  is  defined  to  be  at  the  centre 
of  the  Bragg  cells  and,  correspondingly  at  the  centre  of  their 
images  on  the  detector  array,  the  optical  signals  diffracted  by  the 
Bragg  cells  are  a(t+z/v)  and  b(t-z/v) ,  where  v  is  the  velocity  of 
propagation  of  the  signals  in  the  Bragg  cells.  The  electrical 
signal  produced  by  the  detector  array,  s(t,z),  is  proportional  to 
the  square  of  the  light  distribution 

s(t,z)  =  |a(t+z/v)  +  b(t-z/v)  | 2  (3) 
=  a2(t+z/v)  +  b2(t-z/v)  +  2a(t+z/v)b(t-z/v)  (4) 

A  rigourous  analysis  using  analytical  functions  can  be  found 
in  t32] ♦  Equation  4  is  a  simplification  of  these  results.  This 
response  is  then  integrated  for  a  duration  T  and  the  resulting 
signal  S(T,z)  is 


S(T,z)  =  J  [a2  (t-z/  v)  +  b2  (t-z/v)  +  2a(t+z/v)  b(t-z/v)]dC  (5) 


The  response  of  all  the  individual  elements  of  the  detector 
array  is  called  a  correlogram  (see  Figure  6) .  The  third  term  of 
Equation  3  is  the  correlation  of  the  input  data  a(t)  and  b(t)  . 
This  correlation  term  is  observed  as  a  peak  on  a  pedestal 
(generated  by  the  first  and  second  term  of  Equation  3)  .  The 
presence  of  the  peak  indicates  that  the  two  input  signals  a(t)  and 
b(t)  are  identical,  conversely  the  absence  of  the  correlation  peak 
indicates  that  the  two  inputs  are  different. 

A  photo  of  a  TIC  is  presented  in  Figure  5.  Some  of  the 
advantages  of  this  implementation  is  its  relatively  low  volume  that 
can  be  packaged  into  a  10  inch  high  space  of  a  19  inch  rack.  The 
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Figure  6:  Correlation  peak  from  a  TIC. 
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power  consumption  is  less  than  500  Watts.  The  Bragg  cells  have  a 
30  MHz  bandwidth  centered  at  45  MHz  and  make  the  TIC  well  suited 
for  the  processing  of  spread  spectrum  communications  signals.  The 
detector  array  has  5000  elements  that  can  be  read-out  in  500  /is. 
The  illumination  is  provided  by  a  20  mW  helium-neon  laser. 

The  TIC  is  run  by  a  486  personal  computer  whose  main  functions 
are  to  control  the  read-out  of  the  detector  array  and  to  perform 
the  removal  of  the  pedestal  (Figure  6)  which  facilitates  the 
detection  of  the  peak.  Various  techniques  have  been  proposed  and 
assessed  at  DREO[49]  to  perform  pedestal  removal.  They  all 
involved  the  subtraction  of  correlograms  collected  at  different 
times.  However,  the  time  interval  between  the  collection  of  the 
correlograms  varies  with  the  method.  In  order  to  deal  with  a 
tactical  situation  where  the  signal  environment  is  likely  to  change 
rapidly,  the  method  selected  uses  correlograms  collected  with  the 
minimum  time  interval,  that  is  successive  correlograms.  That 
particular  method  is  described  in[22,35,49] . 

Figure  7  illustrates  the  signal  generation  for  the 
conventional  operation  of  a  TIC  with  the  three  signals  described  in 
section  2.0.  The  input  to  the  TIC  is  made  of  the  addition  of  two 
BPSK  signals  with  chip  rates  of  2MHz  and  5  MHz  respectively,  and 
AWGN  whose  spectra  are  illustrated  in  Figure  3.  The  two  BPSK 
signals  are  added,  and  the  combined  signal  is  split.  Similarly, 
the  AWGN  is  split  and  a  time  delay  of  1.46  /us  is  then  introduced  in 
one  of  the  AWGN  paths.  The  noise  paths  are  then  added  to  the  BPSK 
signals  and  the  two  resulting  signals  are  applied  to  the  Bragg 
cells  of  the  TIC  for  correlation.  The  geometry  of  the  simulated 
signal  sources  is  illustrated  in  Figure  2. 

The  correlation  is  formed  on  the  detector  array  by  an 
integration  of  duration  T,  calculated  to  provide  sufficient 
processing  gain  to  allow  the  detection  of  the  correlation  peaks. 
At  the  end  of  the  integration  time  T,  the  correlation  is  already 
fully  formed  on  the  detector  array  and  is  clocked  out.  This  leads 
to  the  real-time  implementation  of  the  CCC  for  a  particular  cyclic 
frequency  a.  In  the  system  considered  here,  it  takes  a  total  of 
500  /us  to  read  the  data  from  the  5000  elements  of  the  detector 
array  and  access  the  CCC.  However,  the  delay  of  500  /us  is  not  lost 
time,  because  the  detector  array  is  reset  immediately  after  the 
transfer  of  the  results  to  the  device  responsible  for  the  clocking 
out  and  the  next  integration  period  starts  during  the  clocking  out 
of  the  results  of  the  previous  one. 


4.0  IMPLEMENTATION  OF  THE  CYCLIC  AUTOCORRELATION 
WITH  A  TIME-INTEGRATING  CORRELATOR 

The  implementation  of  the  CCC,  R*ls2  (0»  is  suggested  by  the 
interpretation  of  Equation  1  (presented  in  section  2.1)  where  the 
CCC  is  viewed  as  the  correlation  of  one  received  RF  signal  with  a 
frequency  shifted  version  of  the  RF  signal  received  by  another 
antenna.  Remembering  that  a  TIC  produces  the  correlation  of  all  the 
signals  that  are  applied  to  its  Bragg  cells,  the  CCC  can  be 
produced  by  frequency  shifting  the  carrier  frequencies  of  the  two 
RF  input  signals  by  an  appropriate  amount  before  performing  the 
correlation.  A  correlation  peak  appears  when  the  difference  of  the 
carrier  frequencies  of  the  two  signals  is  equal  to  one  of  the 
cyclic  frequencies  contained  in  the  signals. 

A  schematic  of  the  generation  for  a  signal  made  of  the 
addition  of  two  Binary  Phase  Shift  Keyed (BPSK)  with  a  chip  rate  of 
2  MHz  and  5  MHz  and  AWGN  is  illustrated  in  Figure  8.  The  two  BPSK 
signals  are  added,  and  the  combined  signal  is  split  as  in  Figure  7. 
A  time  delay  of  1.46  [ms  is  then  introduced  in  one  of  the  AWGN 
paths.  The  noise  paths  are  then  added  to  the  BPSK  signals.  The 
two  resulting  signals  are  then  mixed  with  two  different  carrier 
frequencies  (master  and  slave)  and  are  applied  to  the  Bragg  cells 
of  the  TIC  for  correlation.  The  conventional  cross  correlation 
r°is2  (t)  is  produced  when  the  two  signals  have  the  same  carrier 
frequency.  The  CCC,  for  a  cyclic  frequency  a,  is  produced  when  a 
frequency  shift  a  has  been  applied  to  the  carrier  frequency  of  one 
of  the  received  RF  signals.  In  this  system,  changes  in  the  carrier 
frequencies  had  to  be  made  by  hand.  The  TDOA  window  of  the 
correlator  used  in  this  work  was  89  /z s,  and  was  detected  by  a 
5000-element  array.  It  is  possible,  without  using  any  special 
algorithm,  to  locate  the  center  of  a  peak  with  a  precision  of  3 
detector  elements.  The  resulting  resolution  for  the  TDOA  is  50  ns. 

4.1  Bandwidth  Limitations 

The  bandwidth  of  the  Bragg  cells  incorporated  in  the  TIC  is  30 
MHz,  centered  at  a  carrier  frequency  of  45  MHz.  In  order  to 
optimize  the  use  of  this  bandwidth,  the  CCC  can  be  implemented  by 
shifting  the  carrier  frequency  of  the  two  signals  by  a/2  and  -a/2 
respectively  rather  than  shifting  the  carrier  of  only  one  signals 
by  a.  This  way,  both  signals  have  a  carrier  frequency  at  the  same 
distance  from  the  center  frequency  of  the  Bragg  cells  and  the 
bandpass  filtering  resulting  frpm  the  frequency  response  of  the 
Bragg  cells  is  the  same  for  both  signals.  Correlation  still 
appears  in  the  CCC  when  the  difference  between  the  carrier  of  the 
two  signals  is  equal  to  one  of  the  cyclic  frequencies  of  the 
signal. 
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Figure  8:  Generation  of  the  input  signals  for  the  CCC. 
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5 . 0  EXPERIMENTAL  RESULTS 


The  purpose  of  the  proof -of -concept  experiments  was  to 
demonstrate  the  capability  of  the  CCC  to  discriminate  between 
signals  that  could  not  be  discriminated  by  a  spectrum  analyser  and 
also  that  optical  TIC  can  perform  real-time  calculation  of  the  CCC 
of  large  bandwidth  spread  spectrum  signals. 

5.1  First  Experiment 

5.1.1  Simulated  Signal  Environment 

The  elements  of  the  signal  environment  were  selected  to 
simulate  a  tactical  situation  involving  three  spread  spectrum 
signals  and  to  emphasize  the  capability  of  the  CCC  to  discriminate 
between  non-stationary  signals  occupying  the  same  bandwidth  and 
having  the  same  TDOA  and  the  same  carrier  frequency.  The  signal 
geometry  is  illustrated  in  Figure  2  and  the  spectra  of  the  signals 
are  illustrated  in  Figure  3.  The  detailed  description  of  the 
signals  and  of  the  signal  generation  have  already  be  presented  in 
Sections  3.0  and  4.0  respectively. 

The  CCC  for  different  cyclic  frequencies  can  be  produced  by 
mixing  different  carrier  frequencies  to  the  signals  before  they  are 
correlated  by  the  TIC  (see  Figure  8)  .  In  our  experiments,  the 
values  of  these  two  carrier  frequencies  are  adjusted  by  hand.  When 
their  difference  is  equal  to  one  of  the  cyclic  frequencies  present 
in  the  input  signals,  a  correlation  peak  appears  in  the  output  of 
the  CCC.  An  automated  scanning  could  be  performed  under  computer 
control. 


5.1.2  Cyclic  Cross  Correlation  Results 

Figure  9  contains  the  CCC  of  the  three  signals  (described 
in  section  4.0)  for  the  cyclic  frequency  varying  between  a=0  and 
a=5  MHz,  by  steps  of  1  MHz.  The  correlograms  produced  by  the  TIC 
are  presented  for  each  cyclic  frequency  a.  The  horizontal  axis  of 
each  correlogram  is  labelled  in  number  of  detecting  elements  called 
pixels.  This  corresponds  to  the  TDOA  after  calibration.  The 
vertical  axis  (with  a  linear  scale)  is  the  amplitude  of  the  CCC. 
Other  interesting  features  can  be  observed  in  these  correlograms: 

1)  Figure  9  a),  for  a=0,  is  the  conventional  cross 

correlation  of  the  input  signals.  In  this  plane,  the 
AWGN  (Noise)  and  the  two  BPSK  (peak  A)  are  separated 
by  their  different  TDOA.  However,  it  is  impossible  to 
separate,  for  a=0,  the  two  BPSK  signals  that  share  the 
same  TDOA  of  0  fxs  and  produce  superposed  correlation 
peaks. 
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Figure  9: 


Output  of  the  optical  correlator  showing 
the  CCC  of  the  sum  of  the  three  signals  of 
Figure  6  for  different  cyclic  frequencies  a 


2)  Figure  9  b)  and  9  d)  illustrates  the  CCC  for  a=l  MHz 
and  a=3  MHz  respectively.  For  these  two  cyclic 
frequencies,  the  AWGN  and  the  BPSK  signals  do  not 
produce  any  correlation  peak. 

3)  For  a=2  the  CCC  contains  one  peak  (Figure  9c) .  Peak 
B  is  associated  with  the  2  MHz  BPSK  signal  and  is 
located  at  t= 0  ns.  An  harmonic  of  the  CCC  for  the  2 
MHz  BPSK  (peak  C)  can  be  found  at  a=4  MHz  (Figure  9e)  . 

4)  Finally,  for  a=5  MHz  (Figure  9f ) ,  peak  D  is  produced 
by  the  5  MHz  BPSK  signal  at  r=0  jus. 

Figure  1  is  a  synthesis  of  the  results  of  Figure  9  and  has 
been  constructed  from  the  experimental  data  collected  in  the  work 
just  described.  It  demonstrates  clearly  that  the  CCC  has  the 
capability  to  separate  BPSK  signals  with  overlapping  spectra, 
sharing  the  same  carrier  frequency  and  the  same  TDOA,  in  the 
presence  of  a  considerable  amount  of  noise  (equal  to  the  power  of 
the  BPSK  signals) . 

5 . 2  Second  Experiment 

5.2.1  Simulated  Signal  Environment 

The  capabilities  of  the  TIC  to  produce  in  real-time  the  CCC 
of  larger  bandwidth  signals  has  been,  explored  by  simulating  the 
signal  environment  illustrated  in  Figure  10.  In  this  case  a  10  MHz 
BPSK,  a  20  MHz  BPSK  and  a  source  of  AWGN  share  the  same  bandwidth, 
carrier  frequency  and  TDOA.  The  spectra  of  the  signals  are 
illustrated  in  Figure  11.  The  schematic  of  the  signal  generation 
is  similar  to  Figure  8  except  that  the  long  cable  producing  the 
time  delay  has  been  removed. 

5.2.2  Cyclic  Cross-Correlation  Results 

The  results  of  the  CCC  for  the  cyclic  frequency  0,  10  MHz, 
20  MHz  and  1  MHz  are  illustrated  in  Figure  12. 

1)  Figure  12  a)  shows  the  CCC  for  a=0.  The  peak  A  is  the 
superposition  of  the  autocorrelation  peak  of  the  three 
signals  present  in  the  environment.  The  peak  B  is  the 
superposition  of  the  autocorrelation  peaks  of  the 
emitters  at  10  and  20  MHz  while  peak  C  is  the 
autocorrelation  peak  of  the  emitter  at  20  MHz. 

2)  Figure  12  b)  shows  the  CCC  for  a=10  MHz.  Peak  D 

indicates  the  presence  of  an  emitter  at  10  MHz. 

3)  Figure  12  c)  shows  the  CCC  for  a=20  MHz.  Peak  E 

indicate  the  presence  of  an  emitter  at  20  MHz. 


18 


Hr.f-  .c  dBm  ATTF-M  lO  dB 

'  1  !  ‘  I  '  '  I 


SPAN 

ao .  ao  mH:: 


center  43  oa  wnr 

RES  BW  l  MHZ 


I  H* r  a  dBm  AT TEH  ID  dB 

,  '  ’  '  l  ' 


VBW  300  kHz 


spa* ;  bo. do 
SWM  P.O  msec 


a)  emitter  1:  10  MHz  BPSK 


b)  emitter  2:  20  MHz  BPSK 


Rf-'P  .0  dBm  ATTEN  *0  dB  +  20 

.  :  :  i  •  ! 


nEF  .0  dBm  ATTEN  10  dB  ♦ 

r  I  *  I 

8/  ; 


span 

80.00  MHZ 


SPAM 

80.00  MHz 


:M  43.00  MHz 

RES  BW  1  mhz 


V0W  300  kHz 


span  nn  on  miiz 

SWM  an  mnnn 


;R  43 . 00  MM y 

RES  BW  i  mhz 


VBW  300  kHZ 


SPAM  RC1  00  MHz 
SWM  20  m*«e 


c)  emitter  3:  AWGN 


d)  sum  of  the  three  emitters. 


Figure  11:  Spectra  of  the  signals  of  the  three  emitters: 
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4) 


Figure  12  d,  shows  the  CCC  for  a  =  11  MHz.  For  this 
cyclic  frequency,  no  peak  is  produced.  Similar  results 
are  obtained  for  all  cyclic  frequencies  that  are 
different  from  10  MHz  and  20  MHz.  It  should  be 
remembered  that  the  bandwidth  of  the  Bragg  cells  used 
in  these  experiments  is  only  30  MHz. 

The  experimental  results  of  Figure  12  are  presented  as  a  3- 
dimensional  computer  generated  picture  on  Figure  13.  This  way  to 
present  the  data  emphasizes  the  capability  of  the  CCC  to  detect  and 
characterize  large  bandwidth  spread  spectrum  signals  even  when  a 
substantial  amount  of  AWGN  is  located  on  the  same  TDOA. 

All  the  experimental  results  presented  in  Figure  9  and  12  were 
obtained  by  letting  the  input  signals  flow  through  the  Bragg  cells 
for  an  integration  time  of  1000  jus.  The  CCCs  for  the  particular 
cyclic  frequency  a  considered  were  fully  formed  on  the  detector 
array  at  the  end  of  the  integration  period  leading  to  the  real-time 
implementation  of  the  CCC.  The  results  were  available  500  /ns  later 
after  being  clocked  out  of  the  detector  array.  The  correlation 
peaks  identified  in  Figure  9  and  12  are  easily  discernible  from  the 
noise  and  from  each  other.  It  would  be  possible  to  calculate  the 
cyclic  spectrum  described  by  equation  2  by  performing  the  Fourier 
transform  of  the  CCC  calculated  by  the  TIC  for  the  desired  value  of 
a. 


5.3  Width  of  the  Correlation  Peaks  on  the 

Cyclic  Frequency  axis  a 

The  width  of  the  CCC  peaks  was  experimentally  measured  by 
recording  the  amplitude  of  the  peak  as  a  function  of  the  cyclic 
frequency.  The  results  are  illustrated  in  Figure  14  for  a  BPSK 
with  a  chip  rate  of  10  MHz.  The  amplitude  of  the  peak  decreases  as 
the  value  of  alpha  gets  away  from  the  cyclic  frequency  of  10  MHz. 
The  amplitude  of  the  peak  becomes  null  500  Hz  away  from  10  MHz. 
However,  a  second  lobe  is  observed  between  500  Hz  and  1000  Hz.  No 
correlation  peak  was  observed  beyond  1000  Hz  away  from  the  chip 
rate  of  the  BPSK.  The  width  of  the  peaks  produced  by  the  other 
chip  rate  used  in  these  experiments  was  identical. 
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Figure  14:  Width  of  the  correlation  peak  on  the  cyclic  frequency 
axis  a,  for  a  BPSK  with  a  10  MHz  chip  rate. 


6 . 0  CONCLUSION 


The  experimental  results  presented  in  this  Technical  Note 
clearly  demonstrate  that  the  CCC  is  a  powerful  algorithm  capable  of 
sorting  spread  spectrum  signals  sharing  the  same  bandwidth,  TDOA 
and  carrier  frequency.  The  characterization  of  the  signals 
performed  by  the  CCC  can  be  done  even  in  presence  of  a  considerable 
amount  of  noise.  It  also  has  been  demonstrated  that  an 
optoelectronic  correlator  is  capable  of  real-time  implementation  of 
the  CCC  for  a  particular  cyclic  frequency.  This  combination  of 
features  makes  an  optoelectronic  implementation  of  the  CCC  a  very 
interesting  candidate  for  the  ESM  of  spread  spectrum  LPI  signals 
processing. 

However,  many  important  issues  have  to  be  studied  and  resolved 
before  attempting  to  design  a  tactical  system.  The  effects  of 
interference  of  various  nature  and  of  multipath  signals  has  to  be 
established.  The  effects  of  non-linearity  in  the  receivers  have  to 
be  assessed.  The  resolution  achievable  on  the  TDOA  axis  r  and  on 
the  cyclic  frequency  a  has  to  be  measured  experimentally  and 
compared  to  operational  requirements.  Finally,  the  dynamic  range 
of  the  system  has  to  be  measured.  It  may  be  advantageous  to 
perform  some  of  these  studies,  at  least  at  the  preliminary  stage, 
by  simulation. 
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